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ABSTRACT 

We study relations between global characteristics of low-redshift (0 < 2 < 1) compact 
star-forming galaxies, including absolute optical magnitudes, H/? emission-line lumi¬ 
nosities (or equivalently star-formation rates), stellar masses, and oxygen abundances. 
The sample consists of 5182 galaxies with high-excitation H II regions selected from the 
SDSS DR7 and SDSS/BOSS DR10 surveys adopting a criterion [O m]A4959/H/3 ^ 1. 
These data were combined with the corresponding data for high-redshift (2 < 2 < 3) 
star-forming galaxies. We find that in all diagrams low -2 and high -2 star-forming 
galaxies are closely related indicating a very weak dependence of metallicity on stellar 
mass, redshift, and star-formation rate. This finding argues in favour of the universal 
character of the global relations for compact star-forming galaxies with high-excitation 
H II regions over redshifts 0 < 2 < 3. 

Key words: galaxies: fundamental parameters - galaxies: dwarf - galaxies: starburst 
- galaxies: ISM - galaxies: abundances. 


1 INTRODUCTION 

Luminosities, stellar masses, and metallicities are among the 
most important properties of galaxies, which are needed for 
a better understanding of their formation and evolution. The 
H /3 emission-line luminosity allows for a direct measure of 
the rate of gas transformation into stars, while stellar mass 
and metallicity are indicators of the galaxy’s past evolution. 
Thus, relations between these parameters and their varia¬ 
tions with redshift may help to analyse the origin and evo¬ 
lution of galaxy populations on cosmological time scales. 

In this respect, the study of star-forming galaxies plays 
a particularly important role. The brightness of these galax¬ 
ies is enhanced by the ongoing star-format ion, allowing us 
to observe them at larger distances. Furthermore, stronger 
emission lines in spectra of these galaxies better constrain 
their metallicity, one of the key galaxy global parameters. 

The Sloan Digital Sky Survey (SDSS) imaging and spec¬ 
troscopy (York et al. 2000) opened a unique opportunity to 
produce the luminosity-metallicity and mass-metallicity re¬ 
lations for large samples of low -2 normal and star-forming 
galaxies. Tremonti et al. (2004) obtained a mass-metallicity 
relation for ~53000 SDSS star-forming galaxies at 2 ~ 0.1. 


They found that the relation is relatively steep for low-mass 
galaxies with stellar masses M* < 10 10 ' 5 Mq, but flattens 
for higher-mass galaxies. The mass-metallicity relation by 
Tremonti et al. (2004) is often used as a benchmark for 
comparisons with similar relations obtained for low -2 and 
high -2 galaxies (e.g. Lee et al. 2006; Maiolino et al. 2008; 
Amorm et al. 2010; Manucci et al. 2010; Maier et al. 2014; 
Steidel et al. 2014; Zahid et al. 2011, 2012, 2013, 2014a,b). 
However, we note that Tremonti et al. (2004) selected all 
star-forming galaxies, including galaxies with low-excitation 
H 11 regions, while only galaxies with high-excitation H II 
regions are present in high -2 samples. They also used cali¬ 
brations by Chariot & Longhetti (2001) for the metallicity 
determination, which were not used for high -2 galaxies. This 
implies that differences between low -2 and high -2 galaxies 
can be introduced, at least in part, by different selection 
criteria and different methods used for the determination 
of global galaxy parameters such as metallicities and stellar 
masses. 

Maiolino et al. (2008) found a strong evolution of the 
stellar mass-metallicity relation with redshift in the sense 
that high -2 (2 ~ 3.5) galaxies are more metal-poor com- 
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pared to low -z galaxies for the same stellar mass. This effect 
is more pronounced for low-mass galaxies. Manucci et al. 
(2010) considered a more general relation between stellar 
mass M*, metallicity, and star-format ion rate (SFR). They 
found that high- 2 ; galaxies are characterised by progres¬ 
sively higher SFRs. Introducing a parameter /x=logM*/M 0 
- cdogSFR(M 0 yr _1 ), Manucci et al. (2010) showed that, in¬ 
dependent of redshift, star-forming galaxies have a uniform 
li - metallicity distribution if a = 0.32. 

On the other hand, there is some evidence that prop¬ 
erties of nearby and high -z star-forming galaxies are similar 
with respect to SFR and metallicity. Heckman et al. (2005) 
have identified nearby (z < 0.3) ultraviolet-luminous galax¬ 
ies (UVLGs) selected from the Galaxy Evolution Explorer 
( GALEX ). These compact UVLGs were eventually called 
Lyman-break analogues (LBAs). They resemble LBGs in 
several respects, implying that mass-metallicity relations for 
low -z and high -z galaxies should be similar. In particular, 
their metallicities are subsolar, and their SFRs of ~ 4 - 25 
M 0 yr -1 are overlapping with those for LBGs. 

Furthermore, Cardamone et al. (2009) selected a sample 
of compact strongly star-forming galaxies from the SDSS, 
which are also similar to LBGs because of their low metal¬ 
licities and high SFRs. Izotov et al. (2011) extracted a sam¬ 
ple of star-forming luminous compact galaxies (LCGs) with 
hydrogen H (3 luminosities L(H/3) ^ 3xl0 40 erg s -1 and H (3 
equivalent widths EW(H/3) ^ 50A from SDSS spectroscopic 
data. These galaxies have properties similar to “green pea” 
galaxies (Cardamone et al. 2009) but are distributed over a 
wider range of redshifts z ~ 0.02 - 0.63. The SFRs of LCGs 
are high at ~ 0.7 - 60 M 0 yr -1 and overlap with those of 
LBGs. 

Izotov et al. (2011, see also Guseva et al. (2009) and 
Zhao et al. (2010)) showed that LBGs, LCGs, luminous 
metal-poor star-forming galaxies (Hoyos et al. 2005), ex¬ 
tremely metal-poor emission-line galaxies at z < 1 (Kakazu 
et al. 2007), and low-redshift blue compact dwarf (BCD) 
galaxies with strong star-formation activity obey a common 
luminosity-metallicity relation. Finally, Izotov et al. (2014a) 
considered luminosity-metallicity and mass-metallicity rela¬ 
tions for SDSS DR7 star-forming galaxies, which are much 
flatter than the Tremonti et al. (2004) relations. They also 
found that low -z galaxies with high SFRs do not deviate 
from the relation established for galaxies with lower SFRs. 

The above discussion implies some inconsistency be¬ 
tween results obtained in different papers and raises the 
question: are luminosity- and mass-metallicity relations be¬ 
tween global parameters of low -z and high -z star-forming 
galaxies different? What is the impact of different se¬ 
lection criteria and methods used for the determination 
of the galaxy metallicities? In this paper we attempt to 
study similarities and differences of the relations for star¬ 
forming galaxies in a redshift range 0 < z < 3. Special 
care has been taken to use strong-line methods to derive 
the metallicity, which were calibrated using the T e -method, 
based on the electron temperature determination from the 
[O iii](A4959+A5007)/A4363 flux ratio (e.g. Izotov et al. 
2011), and thus give consistent gas-phase metallicities. We 
also selected only compact star-forming galaxies with high- 
excitation H 11 regions, for which both methods can be ap¬ 
plied. Furthermore, excitation conditions in those H 11 re¬ 
gions are similar to those in high -z galaxies. 



l°g([NII]A6584/Ha) 

Figure 1. The Baldwin-Phillips-Terlevich (BPT) diagnostic di¬ 
agram (Baldwin et al. 1981). Selected SDSS compact star¬ 
forming galaxies with [O ill] A4959/H/3 ^ 1.0, corresponding to 
log [O iii]A5007/H/3 > 0.5, and detected [N ii]A6584 line are shown 
by blue filled circles. Also plotted are all emission-line galaxies 
from SDSS DR7 (cloud of grey dots). The red solid line from 
Kauffmann et al. (2003) separates star-forming galaxies from ac¬ 
tive galactic nuclei. 

In Sect. 2 we discuss the sample. The technique used 
for the determination of luminosities and stellar masses is 
described in Sect. 3. Different methods of the oxygen abun¬ 
dance determination are discussed in Sect. 4. Relations be¬ 
tween different global parameters for our SDSS sample and 
comparisons with respective relations for high- 2 : galaxies 
with high-excitation H 11 regions are presented in Sect. 5. 
The main results of the paper are described in Sect. 6. 


2 THE SAMPLE 

The sample of compact star-forming galaxies was selected 
from the spectroscopic data base of the SDSS DR7 (Abaza- 
jian et al. 2009) and SDSS/BOSS DR10 (Ahn et al. 2014; 
Dawson et al. 2013; Smee et al. 2013). The criteria were: 
1) only galaxies with angular diameters < 6 " were selected 
to minimize the aperture corrections of the derived lumi¬ 
nosities and stellar masses. As a measure of the angular 
size we used the Petrosian radius R 50 with a 50% flux 
of the galaxy inside it in the SDSS r band. These data 
are available in the SDSS data base. The oxygen abun¬ 
dance is then characterizing the entire galaxy because of 
its compactness; 2) spiral galaxies are excluded. The oxy¬ 
gen abundances in individual H 11 regions are not represen¬ 
tative characteristics for an entire spiral galaxy because of 
the abundance gradients in those galaxies; 3) galaxies with 
emission-line ratios [O m]A4959/H/3 <1.0 were excluded to 
ensure that only galaxies with high-excitation H II regions 
are included, for which the determination of extinction and 
chemical composition is more reliable. We note that most 
of high-redshift emission-line galaxies satisfy the criterion 
[O ill]A4959/H/3 ^ 1.0, and hence have similar excitation 
conditions in their H II regions; 4) galaxies with AGN ac¬ 
tivity are excluded. We used obvious AGN spectral charac¬ 
teristics: broad emission lines (QSOs and Syl galaxies), the 
presence of strong high-ionisation [Ne v] A3426 and He II 
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Figure 2. (a) The 25 // x25 // SDSS image centered on the com¬ 
pact star-forming galaxy J0207+0047. (b) The rest-frame SDSS 
spectrum of J0207+0047. The galaxy is located at a redshift of z 
= 0.54. 



Figure 3. Distribution over redshift z of the SDSS compact star¬ 
forming galaxies from the entire sample. The dotted vertical line 
indicates the average redshift of 0.16. 


A4686 emission lines in high-excitation spectra of Sy2 galax¬ 
ies. Strong high-ionisation nebular [O ill] emission lines are 
simultaneously observed with strong low-ionisation nebular 
[O i], [O n], [N n], and [S n] emission lines in many high- 
excitation Sy2 spectra. LINERs are excluded by a condition 
[O ill] A4959/H/5 < 1.0. The reliability of our selection is 
checked with the diagnostic diagram, proposed by Baldwin 
et al. (1981) (BPT diagram, Fig. 1). All galaxies from our 
sample are located in the region of star-forming galaxies 
(blue dots in Fig. 1). 

Applying these criteria, we selected 5182 galaxies in the 
redshift range 0 < z < 1. The [O ill] A4363 emission line 
was detected at a level higher than lcr in 3607 galaxies from 
our sample, allowing us to obtain oxygen abundances by 
the Te-method. Hereafter this subsample will be called “T e 
sample”. We apply strong-emission-line (SEL) methods to 
derive oxygen abundances for the remaining subsample of 
1575 galaxies. 

The SDSS image of one compact galaxy from our sam¬ 
ple, J0207+0047, is shown in Fig. 2a. This galaxy is located 
at the redshift z = 0.54 and is very compact and is almost 
unresolved. Its angular diameter is only — 1.4 7/ , corre¬ 
sponding to a linear diameter of ~ 30 kpc. However, since 
this galaxy is unresolved, its true linear diameter can be 
much smaller. This is a common property for distant galax¬ 
ies from our sample with z > 0.2 - 0.3. HST images of a 
few compact star-forming galaxies at z ~ 0.1 - 0.35 reveal 
that all of them have linear diameters less than 3 kpc (Car- 
damone et al. 2009; Jaskot &; Oey 2014, Izotov et al. 2015, 
in preparation). The spectrum of J0207+0047 (Fig. 2b) is 
characterised by strong narrow emission lines and blue con¬ 
tinuum. Spectra of all our galaxies have similar character¬ 
istics. They are very different from the spectra of “normal” 
galaxies dominated by numerous absorption lines and much 
redder continua. The spectrum in Fig. 2b is also very differ¬ 
ent from typical spectra of star-forming galaxies with weak 
emission lines, which were selected by Tremonti et al. (2004). 

The redshift distribution of our entire sample with the 
average redshift of 0.16 (dotted vertical line) is shown in 
Fig. 3. 

Spectroscopic data were supplemented with photomet¬ 
ric data in five SDSS bands u , g , r, i, and z for the entire 
galaxy and within the spectroscopic aperture, allowing aper¬ 
ture corrections of the observed fluxes and the determination 
of absolute galaxy magnitudes. 


3 THE DETERMINATION OF GALAXY 
GLOBAL PARAMETERS 

3.1 Parameters of emission lines and dust 
extinction 

We measured emission-line fluxes and equivalent widths us¬ 
ing the IRAF 1 SPLOT routine. The line-flux errors in¬ 
clude statistical errors in addition to errors introduced by 


1 IRAF is the Image Reduction and Analysis Facility distributed 
by the National Optical Astronomy Observatory, which is oper¬ 
ated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science 
Foundation (NSF). 
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the standard-star absolute flux calibration, which we set 
to 1% of the line fluxes. These errors will be later propa¬ 
gated into the calculation of abundance errors. Using the 
observed decrement of several hydrogen Balmer emission 
lines we corrected the line fluxes relative to the H/3 flux 
for two effects: (1) reddening adopting the extinction curve 
of Cardelli et al. (1989) and (2) underlying hydrogen stel¬ 
lar absorption that is derived simultaneously by an iterative 
procedure as described in Izotov et al. (1994). The extinc¬ 
tion coefficients are defined as C(R/3) = 1A7E(B—V), where 
E(B -V) = A(V)/ 3.2 (Aller 1984). 

The spectroscopic data were corrected for the aperture 
using the relation 2.5 r( ' app ' )-r , where r and r(app) are the 
SDSS r-band total magnitude and the magnitude within the 
round spectroscopic aperture, respectively. The diameters of 
spectroscopic apertures are 3" and 2" for the SDSS DR7 
and SDSS/BOSS spectra, respectively. Selected galaxies are 
compact. Therefore the aperture corrections of fluxes are 
small, rarely exceeding the value of 2. 

3.2 The absolute g-band magnitudes and H /3 
luminosities 

The extinction-corrected absolute g-band magnitudes M g 
and extinction- and aperture-corrected Hluminosities 
L(R/3) were obtained respectively from the SDSS extinction- 
corrected g magnitude for the entire galaxy and extinction- 
and aperture-corrected H /3 emission-line flux. The distance 
is derived from the redshift. For distance determination we 
used the relation D — from Refsdal et al. 

(1967), where the Hubble constant H 0 = 67.3 km s -1 Mpc -1 
and cosmological parameters, 9 m = 0.273 and Ua = 0.682, 
were obtained from the Planck mission data (Ade et al. 
2014). The equivalent widths EW(H/3) were correspondingly 
reduced to the rest frame. 


3.3 Spectral energy distribution (SED) and the 
galaxy’s stellar masses 

The stellar mass is one of the most important global galaxy 
characteristics. For its determination we follow the prescrip¬ 
tions described by Guseva et al. (2006, 2007) and Izotov et 
al. (2011, 2014a,b). 

The method is based on fitting a series of model SEDs 
to the observed one and finding the best fit. It consists of 
the following. The fit was performed for each SDSS spec¬ 
trum over the entire observed spectral range of AA3900- 
9200A for SDSS DR7 galaxies and of AA3600-10300A for 
SDSS/BOSS galaxies. As each SED is the sum of both 
stellar and ionised gas emission, its shape depends on the 
relative contribution of these two components. In galaxies 
with high EW(H/3) > 50A, the ionised gas continuum is 
strong and is subtracted before fitting the stellar SED. Izo¬ 
tov et al. (2011) noted the importance of gaseous continuum 
subtraction, otherwise stellar masses of galaxies with high 
EW(H/3)’s in their spectra would be overestimated by a fac¬ 
tor of ~ 3 or more. Hence, for the spectra of galaxies with 
EW(H^) > 50A, SED fitting and subtraction of the ionised 
gas emission in the continuum and emission lines are manda¬ 
tory. Photometric data are not sufficient in these cases. 

We carried out a series of Monte Carlo simulations to 


reproduce the SED of each galaxy in our sample. To derive 
the stellar SED, we use a grid of instantaneous burst SEDs in 
a wide range of ages from 0.5 Myr to 15 Gyr calculated with 
package PEGASE.2 (Fioc & Rocca-Volmerange 1997). We 
adopted a stellar initial mass function with a Salpeter slope, 
an upper mass limit of 100 Mq, and a lower mass limit of 
0.1 Mq. Then the SED with any star-formation history can 
be obtained by integrating the instantaneous burst SEDs 
over time with a specified time-varying star-formation rate. 

The SED of the gaseous continuum was taken from Aller 
(1984). It included hydrogen and helium free-bound, free- 
free, and two-photon emission. 

The star-formation history in each galaxy is approxi¬ 
mated assuming a recent short burst with age t y < 10 Myr, 
which accounts for the young stellar population, and a prior 
continuous star formation for the older stars during the time 
interval between ti and tf (£/ < U and zero age is now). The 
contribution of each stellar population to the SED was pa¬ 
rameterized by the ratio b—M y /M 0 , where M y and M 0 are 
respectively the masses of the young and old stellar popula¬ 
tions. Then the total stellar mass is M*= M y + M 0 . 

For each galaxy, we calculated 10 4 Monte Carlo models 
by randomly varying t y ,U,tf, and b. The best modelled SED 
was found from y 2 minimization of the deviation between 
the modelled and the observed continuum in five wavelength 
ranges, which are free of the emission lines and residuals of 
the night-sky lines. Typical stellar mass uncertainties for our 
sample galaxies are of ~ 0.1 - 0.2 dex. 


4 OXYGEN ABUNDANCES 

One of the most reliable methods of the oxygen abundance 
determination is a T e -method, which is based on the electron 
temperature derived from the ([O ill] A4959+A5007)/A4363 
flux ratio. This method provides a good measure of the 
metallicity below about 12+logO/H < 8.3 (e.g. Pilyugin 
2001; Pettini & Pagel 2004; Stasinska 2005). However, at 
super-solar metallicities the T e -method tends to underesti¬ 
mate significantly the true metallicity (Stasinska 2005). 

We apply the T e -method for the T e sample of the galax¬ 
ies with [O ill] A4363 emission-line intensity measured at the 
level better than 1 a. As it can be seen below, all these galax¬ 
ies satisfy condition 12+logO/H < 8.3. To determine oxygen 
abundances, we generally followed the procedures of Izotov 
et al. (1994, 1997, 2006), which are discussed in detail by 
e.g. Izotov et al. (2011, 2014a,b). 

Recently, Palay et al. (2012) obtained new effective col¬ 
lision strengths for the [O ill] transitions, which in the tem¬ 
perature range 10000 - 20000K are higher by ~20% for the 
transitions [O ill]A4959, A5007, but are higher by only a few 
percent for the transition [O ill]A4363. These updates in 
the atomic data would result in lower electron temperatures 
by ~500 - 700K and hence in higher oxygen abundances 
(Nicholls et al. 2013) as compared to those obtained with 
the Izotov et al. (2006) iterative procedure. However, this 
effect will increase 12 + logO/H by not more than ~ 0.05 
dex. 

Nicholls et al. (2013) and Dopita et al. (2013) consid¬ 
ered the possibility that electrons in H II regions may not be 
in thermal equilibrium, i.e. their velocity distribution devi¬ 
ates from a Maxwellian distribution. The overall effect will 
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Figure 4. A comparison of oxygen abundances 12 + logO/H, derived by the T e -method (this paper) and different strong-line methods 
from Chariot k Longhetti (2001) for our T e sample. In each panel, the emission-line ratios used by strong-line methods are given in 
parentheses. Galaxies with O 32 =[0 ill]A500T/[O ll]A3727 > 4.5 and < 4.5 are given by magenta filled circles and cyan open circles, 
respectively. Solid lines show the one-to-one relations. 




Figure 5. A comparison of oxygen abundances 12 + logO/H derived by the T e -method (this paper) and strong-emission-line (SEL) 
methods proposed by Pettini k Pagel (2004), which utilize the (a) [N n]A6584/H<a and (b) [O iii]A5007/[N n]A6584 line-intensity ratios. 
Shown are the same galaxies by the same symbols as in Fig. 4. Solid lines denote the one-to-one correlations. 


be an increase of the oxygen abundance. However, for low- 
metallicity high-excitation H II regions selected in this study 
the increase of 12 + logO/H is unlikely to be higher than 0.2 
dex (Dopita et al. 2013). We note that a deviation from the 
Maxwellian distribution is found in the solar wind, but not 
in H 11 regions. Therefore, we use equations by Izotov et al. 
(2006), which assume a Maxwellian distribution of electrons. 

The SDSS spectra of most star-forming galaxies are 
noisy. Then the [O ill] A4363 line is not detected or is mea¬ 
sured with insufficient accuracy. In these cases, strong-line 
methods are commonly used to derive oxygen abundances 
instead of the T e -method. There are two different approaches 
to obtain calibrations for the metallicity by using strong-line 
intensities. One approach is based on deriving the line inten¬ 
sities by producing a grid of H 11 region models with varying 
metallicity, ionisation parameter, and other input parame¬ 
ters (e.g. Chariot k Longhetti 2001; Kewley k Dopita 2002; 
Kobulnicky k Kewley 2004; Dawson et al. 2013). The sec¬ 
ond approach uses the empirical relations between certain 
combinations of strong-line intensities and metallicity, which 


is derived by the T e -method (e.g. Alloin et al. 1979; Pagel 
et al. 1979; Pilyugin 2001; Pettini k Pagel 2004; Pilyugin k 
Thuan 2005; Stasinska 2006; Pilyugin et al. 2010). 

The advantage of all these SEL methods is that they are 
based on the strongest emission lines and thus can be ap¬ 
plied to large galaxy samples. The disadvantage is that most 
of these methods depend not only on the metallicity, but 
also on the ionisation parameter. The detailed comparison 
of strong-line methods (e.g. Kewley k Dopita 2002; Kewley 
k Ellison 2008; Maiolino et al. 2008) shows that they are 
not consistent, non-monotonic (e.g. the 7^23-method, Pagel 
et al. 1979; Kewley k Dopita 2002), depend on parameters 
different from metallicity and result in oxygen abundances, 
which can deviate by ~ 0.5 dex or more. 

To illustrate the problem we compare the oxygen abun¬ 
dances derived for the T e sample by the T e -method and 
different strong-line calibrations by Chariot k Longhetti 
(2001), which are based on population-synthesis and H 11 
region models. Calibrations by Chariot k Longhetti (2001) 
were used by Tremonti et al. (2004) to produce luminosity- 
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Figure 6. (a) - (c) Relations between the oxygen abundances 12+logO/H derived by the T e -method and the line intensity ratios 
[N n]A6584/Ho;, [O iii]A 5007/[N ii] A6584, and [O ill] A500T/[O n]A3727, respectively. Shown are only the galaxies from the T e sample. 
Solid and dotted lines represent maximum-likelihood fits to the data and la dispersions, respectively, dashed lines in (a) and (b) refer to 
the relations obtained by Pettini & Pagel (2004). (d) The relation between the oxygen abundances derived using the 03N2 (Eq. 2) and 
0302 (Eq. 3) calibrations. The solid line is the one-to-one correlation. Green open circles and brown filled circles in all panels correspond 
to galaxies with the redshifts z < 0.3 and z > 0.3, respectively. 


metallicity and mass-metallicity relations for a large sample 
of nearby SDSS star-forming galaxies. The Tremonti et al. 
(2004) relations are often used as templates of the present- 
day luminosity-metallicity and mass-metallicity relations for 
comparison with similar relations at higher redshifts. 

We compare in Fig. 4 oxygen abundances derived by us¬ 
ing the Te-method and four out of seven strong-line calibra¬ 
tions by Chariot & Longhetti (2001), which are based on dif¬ 
ferent emission-line intensities and their ratios. The sample 
is split into two parts with the extinction-corrected ratio O 32 
= [O hi] A5007/[0 11 ]A3727 ^ 4.5 (small magenta filled cir¬ 
cles) and <4.5 (small cyan open circles). The parameter O 32 
is an increasing function of the ionisation parameter at given 
metallicity. Among the Chariot & Longhetti (2001) cali¬ 
brations the largest number of emission lines ([O 11 ] A3727, 
[O hi] A5007, [N 11 ] A6584 and [S 11 ] A6717, 6731) is used in 
the calibration, which Chariot & Longhetti (2001) labelled 
as case A. We do not consider two calibrations by Chariot & 
Longhetti (2001), which use the blend of Ha and [N 11 ] A6584 
emission lines. This is because the spectral resolution of the 
SDSS spectra is sufficient to separate these two lines. We 
also excluded from the consideration one calibration, which 
is based on the discontinuity parameter D 4000 at A ~ 4000A. 
This calibration is suitable for galaxies with absorption lines 
in their spectra. 

It is seen in Fig. 4a that the oxygen abundance at low 
metallicities (12 + logO/H (T e ) < 8.1) derived by the strong¬ 
line method (case A) is by ~ 0.2 - 0.3 dex higher than that 
derived by the T e -method. At higher metallicities the oxy¬ 
gen abundances derived by the two methods are very differ¬ 
ent. The comparison of the T e -method with other strong-line 


methods by Chariot & Longhetti (2001) (Fig. 4b - 4d) also 
indicates large inconsistencies. 

Therefore, we conclude that the oxygen abundances, 
and consequently, the luminosity-metallicity and mass- 
metallicity relations for nearby SDSS star-forming galaxies, 
obtained by Tremonti et al. (2004), may not be correct, be¬ 
cause they use a set of Chariot & Longhetti (2001) internally 
inconsistent calibrations. 

A more promising way is to use empirical relations be¬ 
tween line intensities and oxygen abundances derived by the 
Te-method in real star-forming galaxies (e.g. Pettini & Pagel 
2004). These relations can be obtained only for galaxies with 
high-excitation H II regions, where the [O ill] A4363 emission 
line can be detected. Our T e sample satisfies these condi¬ 
tions. In Fig. 5 we compare oxygen abundances obtained 
by the T e -method with those, which were obtained by using 
strong-line calibrations by Pettini &; Pagel (2004). It is seen 
that abundances derived by different methods are in agree¬ 
ment despite the larger dispersion and small offset of ~ 0.1 
at 12+logO/H > 8.0. 

Using only SDSS data for the large T e sample we mod¬ 
ified the Pettini & Pagel (2004) calibrations (Fig. 6) and 
obtained the maximum-likelihood linear relations 

12 + log2 = 8.64 ± 0.02 + (0.47 ± 0.02) x N2, (1) 

H 

12 + log2 = 8.59 ± 0.02 - (0.28 ± 0.01) x 03N2, (2) 

H 

where N2 = log([N ii]A6584 / Ha), OSN2 = log([0 iii]A5007 
/ H J3) - log([N n]A6584 / Ha). 

It is seen in Fig. 6 that relations by Pettini & Pagel 
(2004) (dashed lines) obtained for a smaller sample of H II 
regions, but for a larger range of oxygen abundances with 
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Figure 7. The relations between extinction-corrected absolute 
SDSS g-band magnitude M g and oxygen abundance 12+logO/H. 
The galaxies with the [O ill] A4363 emission-line detection at a 
level higher than 4a shown by brown filled circles and the rest 
of the sample by cyan open circles. The solid line provides a lin¬ 
ear maximum-likelihood fit to the data and dotted lines show la 
dispersions. The dashed line is a linear maximum-likelihood fit 
to the data shown by brown symbols. The steep dash-dotted line 
indicates the relation by Tremonti et al. (2004). The location of 
Lyman-break galaxies (LBGs) at z ~ 3 by Pettini et al. (2001) is 
marked by a grey rectangle. 


inclusion of high-metallicity H II regions, are not very dif¬ 
ferent from the relations Eqs. 1 and 2 (solid lines) obtained 
in this paper. However, the difference is smaller (^0.1 for 12 
+ logO/H) in the case of OSN2 calibrations (Fig. 6b). The 
galaxies from our entire sample were selected using the same 
selection criteria, as the T e sample, i.e. H n regions in these 
galaxies are of high excitation. The only difference is that 
the [O ill] A4363 was measured with insufficient accuracy or 
not detected in many galaxies. Therefore, calibrations Eqs. 
1 and 2 can be applied for the entire sample. 

However, there are 638 galaxies from the entire sam¬ 
ple, including 264 galaxies from the T e sample, where the 
[N n] A6584 emission line has not been measured because 
of its weakness or because this line is outside the SDSS 
wavelength range in high-redshift galaxies with z > 0.40 
of the SDSS DR7 survey and * > 0.55 of the SDSS/BOSS 
survey. In all these cases, with exception of the galaxies 
from the T e sample, we use the relation between the oxy¬ 
gen abundance 12+logO/H and the parameter 0302 = 
log([0 iii]A 5007/[0 n] A3727) derived for the T e sample (Fig. 
6c) 

12 + log2 = 8.18 ± 0.01 - (0.43 ± 0.02) x 03 02. (3) 

H 

We compare in Fig. 6d oxygen abundances derived from 
Eqs. 2 and 3. It is seen that 12 +logO/H obtained with 
the both 03N2 and 0302 calibrations are nearly identical. 
Furthermore, in all panels of Fig. 6 there is no evident offset 
between the galaxies with low redshifts z ^ 0.3 (green open 
circles) and higher redshifts z > 0.3 (brown filled circles), 
implying that all three relations Eqs. 1-3 can be applied 
in the entire range of redshifts. 

Summarising, we derive oxygen abundances (1) by the 
Te-method in 3607 galaxies, (2) by using Eq. 3 in 638 - 264 
= 374 galaxies, and (3) by using Eq. 2 in the remaining 
galaxies. 


5 RESULTS 


5.1 Absolute g-band magnitude - metallicity 
relation 

The extinction-corrected absolute g-band magnitude - 
metallicity relation for the entire sample of star-forming 
galaxies is shown in Fig. 7. The maximum-likelihood linear 
regression, 


12 + log 2 = -(0.041 ± 0.002)M 9 + 7.27 ± 0.03, 

H 


(4) 


is shown by a solid line and la deviations by dotted lines. 
For comparison, we show by dashed line the maximum- 
likelihood linear regression for the T e subsample with the 
[O ill] A4363 emission line detected at the level higher than 
4cr. The difference between these two relations is very small, 
less than 0.05 in 12 + logO/H. This indicates that decreasing 
the detection limit for the [O ill] A4363 emission line to la 
and including galaxies with 12 + logO/H derived by 03N2 
or 0302 calibrations leads to a higher dispersion by a fac¬ 
tor of ~2 but does not introduce significant offsets in 12 + 
logO/H. 

The relation in Eq. 4 is slightly flatter than that ob¬ 
tained by Guseva et al. (2009) and Izotov et al. (2011). 
However, we note that Guseva et al. (2009) and Izotov et al. 
(2011) considered not only SDSS galaxies but also galaxies 
from other samples with a large fraction of low-metallicity 
galaxies. Our sample is more uniform and includes only 
SDSS compact galaxies. The M g - 12+logO/H relation for 
this sample is very similar to that derived by Izotov et al. 
(2014a) for a much smaller sample of SDSS compact star¬ 
forming galaxies. 

Our g-band absolute magnitude - metallicity relation is 
much flatter than that obtained by Tremonti et al. (2004) for 
a sample of ~ 53000 star-forming galaxies selected from the 
SDSS (dash-dotted line in Fig. 7). This difference arises pri¬ 
marily because of different methods applied for the metallic¬ 
ity determination and because of the selection criteria used 
in both studies. Our sample consists of compact galaxies 
with high-excitation H n regions, and we applied strong-line 
methods, which give oxygen abundances consistent with the 
Te-method, at variance with Chariot &; Longhetti (2001) 
calibrations used by Tremonti et al. (2004). Furthermore, 
most of the star-forming galaxies from the Tremonti et al. 
(2004) sample are galaxies with low-excitation H II regions, 
for which the determination of metallicities is more uncer¬ 
tain compared to galaxies with high-excitation H n regions. 
Since empirical calibrations inferred from observations are 
not available for low-excitation H II regions because of the 
very weak and undetectable [O ill] A4363 emission lines, we 
do not consider these galaxies. 

On the other hand, the location of star-forming galaxies 
at 0 < z < 1 with high luminosities is similar to that for the 
Ly-break galaxies (LBGs) by Pettini et al. (2001), shown by 
a box in Fig. 7. This is consistent with a universal character 
of the M g - 12+logO/H relation for star-forming galaxies 
with high-excitation H n regions in the entire redshift range 
0 < z < 3. 
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Figure 8. The H (3 luminosity-metallicity relation. The upper axis indicates the scale of the star-formation rate, calculated with Kennicutt 
(1998) calibration. Galaxies with the equivalent widths EW(H^) <100A and ^100A are shown by the open cyan circles and filled magenta 
circles, respectively. The H/3 luminosities in panel (a) are corrected for extinction and aperture, while they are additionally reduced to 
zero starburst age in panel (b). Solid lines in both panels are maximum-likelihood linear fits to all data and dotted lines are lcr dispersions. 
The dashed lines are maximum-likelihood linear fits only for the galaxies with EW(H/3) ^100A shown by magenta symbols. The location 
of star-forming galaxies with redshifts 2.0 < z < 2.6 by Steidel et al. (2014) and with redshifts z > 2 by Cullen et al. (2014) is shown by 
grey rectangles. Star-forming galaxies with redshifts 2 : ~ 3.4 by Troncoso et al. (2014) are indicated by dark asterisks. 
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Figure 9. The stellar mass - metallicity relations for (a) - (b) the T e and (c) - (d) the entire sample. Brown symbols in (a) correspond 
to the galaxies, where [O ill] A4363 emission is detected at a level better than 4cr, while cyan symbols are for the rest of the galaxies. 
Galaxies shown with magenta and green symbols in (b) - (d) have SFR ^ 10 M@ yr -1 and < 10 Mq yr -1 , respectively. Thick solid 
lines in (a) - (c) represent maximum-likelihood linear fits to the data in respective panels and dotted lines indicate lcr dispersions. The 
thick dashed line in (c) is the maximum-likelihood quadratic fit to the data. The thick solid line in (d) is the maximum-likelihood linear 
fit to the data in (a). The relation by Tremonti et al. (2004) is shown by thin dashed lines and dot-dashed lines are lcr dispersions. The 
location of stacked data for star-forming galaxies at 2 : > 2 by Cullen et al. (2014) and of data for ninety six star-forming galaxies at 
2.0 < 2 ; < 2.6 by Steidel et al. (2014) are shown by dark-grey filled stars and open circles, respectively. Star-forming galaxies at 2 ~ 
3.5 by Maiolino et al. (2008) and at 2 ; ~ 2 by Maier et al. (2014) are shown by dark-grey filled squares and crosses, respectively. Thirty 
star-forming galaxies with redshifts 2 : ~ 3.4 by Troncoso et al. (2014) are indicated by asterisks. The relation for 2 : ~ 2.3 by Sanders et 
al. (2015) is shown by thick yellow solid line. Oxygen abundances for all high- 2 : galaxies are calculated with 03N2 or 0302 calibrations 
from this paper (Eqs. 2-3). (d) Same as in (c), but oxygen abundances are derived using the Chariot & Longhetti (2001) strong-line 
method (case A). 


5.2 H (3 luminosity - metallicity relation 

H (3 emission is a signature of the youngest stellar population 
in star-forming galaxies and the H/3 luminosity allows us 
to estimate the SFR using, e.g., the relation by Kennicutt 
(1998) between SFR and Ha luminosity. However, we note 
that the use of the relation by Kennicutt (1998) is not well 


justified for starburst galaxies because of uncertainties in the 
adopted duration of the burst. The H/3 and Ha luminosities 
are more reliable characteristics, because they are directly 
derived from observations. 

The UV, H (3 and mid-infrared luminosities of compact 
star-forming galaxies are tightly correlated (e.g. Izotov et 
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al. 2014a) indicating that the radiation of young stellar pop¬ 
ulations is a dominant source of dust heating and of galaxy 
emission in the UV range and H (3 line. Consequently, SFRs 
derived by different methods for star-forming galaxies in a 
wide range of redshifts can be compared. 

In Fig. 8a we show the extinction- and aperture- 
corrected H [3 luminosity - oxygen abundance relation for our 
entire sample. The upper axis denotes the SFR calculated 
with the Kennicutt (1998) relation. The box indicates the lo¬ 
cation of star-forming galaxies with redshifts 2.0 < z < 2.6 
by Steidel et al. (2014) and with z > 2 by Cullen et al. 
(2014). Star-forming galaxies with z ~ 3.4 by Troncoso et 
al. (2014) are indicated by asterisks. Strong emission lines 
with [O III] A4959/H/3 > 1 are present in spectra of most 
galaxies from their samples. Therefore, strong-line methods 
calibrated with the T e -method can be applied. Steidel et al. 
(2014) compared SFRs obtained for their sample from the 
UV continuum luminosities and Ha luminosities and found 
good agreement. They derived oxygen abundances by using 
strong-line 03N2 calibration, similar to that applied for our 
sample. To derive metallicities Cullen et al. (2014) and Tron¬ 
coso et al. (2014) used empirical calibrations by Maiolino et 
al. (2008), which also were calibrated with the use of galaxy 
samples, where oxygen abundances were derived by the Te- 
method. However, for consistency reasons we selected only 
high-z galaxies with [O hi] A4959/H /3 ^ 1 and recalculated 
their oxygen abundances using Eqs. 2 or 3. 

Despite the noisy data for high-z galaxies there is a good 
coincidence in Fig. 8a between the luminous tail in the dis¬ 
tributions of our galaxies and the high -z galaxies, indicating 
the universal character of the L(H/3) - 12+logO/H relation 
for star-forming galaxies with high-excitation H n regions. 

However, we note the offset between the distributions 
of the galaxies with EW(H/3) ^ 100A (magenta symbols and 
dashed line) and of the galaxies with EW(H/3) < 100A (cyan 
symbols and solid line). This difference is primarily due to 
the fact that starbursts with high EW(H/3) are on average 
younger. The H/3 emission is powered by the ionising ra¬ 
diation and its luminosity strongly declines with the age of 
starburst. According to the Starburst99 models (Leitherer et 
al. 1999) the number of ionising photons and, respectively, 
the H/3 luminosity produced by a burst is constant during 
the first 3 Myr. After that time it declines and is lower by 
a factor of ~ 50 at a starburst age of 10 Myr. On the other 
hand, the luminosity of non-ionising radiation at ~ 1500 A 
during the same time interval is decreased only by a factor 
of ~ 2. Since galaxies are observed with different starburst 
ages, their H (3 luminosities, for consistency, should be re¬ 
duced to zero age. The H /3 equivalent width can be used for 
this, because it depends on the starburst age. We use Star- 
burst99 models to derive the correction of the H/3 luminosity 
for the starburst age, 

A log L(H/3) = 2.553 - log[EW(H/3)] (5) 

for logEW(H (3) < 2.553, otherwise A log L(H/3) = 0. The 
relation Eq. 5 is correct for instantaneous bursts and may 
overestimate the correction of the H [3 luminosity for the star- 
burst age if underlying older stellar population is present. 
Our SED fitting shows that the contribution of an older stel¬ 
lar population to the stellar continuum near the H/3 emission 
line in galaxies with strong emission lines varies in a wide 
range of ~ 5% - 50%. This results in an overestimation of 


logL(H^) by < 0.3. The effect is stronger for galaxies with 
low EW(H/3). 

In Fig. 8b we show the relation between the H/3 luminos¬ 
ity To (H/3) reduced to zero age and the oxygen abundance. 
It is seen that the differences between maximum-likelihood 
linear fits for our two samples (solid and dashed lines) are re¬ 
duced and distributions of data are in better agreement com¬ 
pared to Fig. 8a, implying that the age correction should be 
applied for consistent comparison of H/3 and Ha luminosities 
and star-format ion rates in star-forming galaxies. 


5.3 Stellar mass - metallicity relation 

In Fig. 9a we show the extinction- and aperture-corrected 
stellar mass - metallicity relations for the T e sample. The 
sample is split into two subsamples with [O ill] A4363 emis¬ 
sion line detected at the level better than 4a (small filled 
brown circles) and in the range la - 4a (small open cyan 
circles). The data for the T e sample can be fitted by a linear 
relation, as shown in Fig. 9a by the solid line, 

12 + log 2 = (0.107 ± 0.004) log + 7.10 ± 0.04. (6) 

H Mq 

The dispersion for the subsample with weaker 
[O ill] A4363 lines is characterised by a factor of ~2 higher 
value, but no offset relative to the subsample with stronger 
lines is present. Similarly, the subsample with high star- 
formation rate SFR b 10 M® yr _1 (small filled magenta 
circles in Fig. 9b) follows the same relation as the subsam¬ 
ple with low SFR <10 M® yr -1 (small open green circles), 
indicating that there is no obvious dependence of the metal¬ 
licity on SFR. 

In Fig. 9c the stellar mass - metallicity relation is shown 
for the entire sample. Adding the galaxies with the oxygen 
abundance derived with strong-line methods results in a flat¬ 
tening of the relation at M* ^ 10 10 M®. Linear regression 
yields 

12 + log 2 = (0.089 ± 0.004) log |L + 7.28 ± 0.04. (7) 

H AZ® 

However, better agreement can be obtained with the 
quadratic stellar mass - metallicity relation 

12 + log 2 = log ^0.023 log ^ + 0.466^) + 5.72 (8) 

shown in Fig. 9c by a thick dashed line. These relations are 
similar to the relation obtained by Izotov et al. (2014a) for 
a smaller sample of SDSS galaxies. 

On the other hand, the stellar mass-metallicity relations 
in Fig. 9a - 9c (solid lines) are much flatter and are shifted to 
lower metallicities compared to that obtained by Tremonti 
et al. (2004) (thin dashed lines). In part, this difference may 
be caused by different selection criteria. However, we argue 
that the main cause of the difference is the use of different 
methods for the oxygen abundance determination, as dis¬ 
cussed above. We show in Fig. 9d the stellar mass - metal¬ 
licity relation for the entire sample, same as in Fig. 9c, but 
with oxygen abundances derived by the Chariot V Longhetti 
(2001) calibration, which has been used by Tremonti et al. 
(2004). It is seen that agreement with the Tremonti et al. 
(2004) relation is much better if the Chariot V Longhetti 
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Table 1 . Average oxygen abundances for the entire sample in the four stellar mass intervals 


SFR 


12+logO/H a 


{Mq yr- 1 ) 

9.0 - 9.2 b 

9.4 - 9.6 b 

10.0 - 10.2 b 

10.4 - 10.6 b 

0.1 - 1 

1 - 5 

5 - 10 

10 - 20 

20 - 50 

> 50 

8.10±0.18 (205) 
8.10±0.16 (147) 
8.11±0.15 ( 67) 
8.12±0.16 ( 52) 
8.12±0.15 ( 27) 

8.12±0.16 (142) 
8.14±0.15 (163) 
8.14±0.13 ( 69) 
8.12±0.20 ( 42) 
8.07±0.18 ( 46) 
8.08±0.15 ( 22) 

8.12±0.16 ( 25) 
8.12±0.15 ( 24) 
8.15±0.16 ( 30) 
8.13±0.11 ( 28) 
7.98±0.20 ( 14) 

8.14±0.09 ( 6) 

8.15±0.07 ( 4) 

8.09±0.13 ( 5) 

8.09±0.21 ( 15) 
8.12±0.13 ( 9) 


a The number of galaxies used for averaging is shown in parentheses. 
b interval of log M* /Mq . 


(2001) calibration is used. However, this calibration is not 
consistent with the T e -method, as it was discussed above. 

It was argued in some recent papers (e.g. Manucci et al. 
2010; Lara-Lopez et al. 2010; Hunt et al. 2012; Andrews & 
Martini 2013) that star-forming galaxies with higher SFRs 
are systematically more metal-poor on the stellar mass- 
metallicity diagrams. In particular, Manucci et al. (2010) 
proposed a fundamental metallicity relation (FMR) between 
stellar mass, gas-phase metallicity and star formation rate. 

We checked whether the tendency of lower metallicity 
in galaxies with higher SFR is present for galaxies from our 
sample. To exclude the dependence of the oxygen abundance 
on the stellar mass we calculated average 12 + logO/H for 
galaxies from the entire sample in the four narrow ranges 
of stellar masses and for different ranges of SFR (Table 1). 
We find no evidence that the galaxies with high SFRs are 
systematically more metal-poor, confirming conclusions by 
Izotov et al. (2014a). We also applied the x 2 minimization 
technique for samples shown in Fig. 9b in the plane [fi a — 
log(M*/M q) - a;logSFR(M 0 yr -1 )] - metallicity and find a 
~ 0, i.e. no dependence on SFR is detected. 

Manucci et al. (2010) argued that introducing with 
a = 0.32 eliminates any redshift evolution up to z ~ 2.5 
yielding the same - metallicity relation independent of 
redshift (see also Sect. 1). To verify this conclusion with our 
sample we show in Fig. 9 star-forming galaxies at z > 2 by 
Cullen et al. (2014) (stars), star-forming galaxies at 2.0 < 
z < 2.6 by Steidel et al. (2014) (filled circles), star-forming 
galaxies at redshifts z ~ 3.5 by Maiolino et al. (2008) (filled 
squares), star-forming galaxies with z ~ 2 by Maier et al. 
(2014) (crosses), and star-forming galaxies with z ~ 3.4 by 
Troncoso et al. (2014) (asterisks). The relation for galaxies 
with z ~ 2.3 by Sanders et al. (2015) is shown by thick yel¬ 
low solid line. For homogeneity, only high-z galaxies with 
[O hi] A4959/H/5 ^ 1 and oxygen abundances derived by 
Eqs. 2 or 3 are shown. It is seen that high-z galaxies are 
closely related to lower-z counterparts and do not indicate 
any appreciable dependence of the mass - metallicity rela¬ 
tion on redshift. However, we note that the high-z data by 
Maiolino et al. (2008), Troncoso et al. (2014), and Maier et 
al. (2014) are offset to lower metallicities compared to the 
data by Cullen et al. (2014) and Steidel et al. (2014). This 
offset may indicate some evolution with the redshift, because 
Maiolino et al. (2008) and Troncoso et al. (2014) considered 
galaxies at higher redshifts (but not Maier et al. 2014). The 
mass-metallicity relation for high-z galaxies is considerably 
flatter and is in better agreement with the relation for low-z 
galaxies in the mass range log M*/Mq — 8.5 - 11, if only 


data by Cullen et al. (2014) and Steidel et al. (2014) are 
considered. 

5.4 H (3 luminosity - stellar mass relation 

Finally, we consider the relation between the H/3 luminosity 
(or, equivalently, the SFR) and stellar mass, both corrected 
for extinction and aperture. At variance with relations dis¬ 
cussed above this relation does not depend on the technique 
used for the determination of metallicity. Fig. 10a shows the 
relation between the extinction- and aperture-corrected H (3 
luminosities (or SFRs) and stellar masses for the entire sam¬ 
ple. The linear maximum-likelihood fit to the data is shown 
by a solid line. The sample is characterised by very high spe¬ 
cific star formation rates sSFR = SFR/M* ranging from ~ 
0.1 Gyr -1 to ~ 100 Gyr -1 (dotted lines), which are among 
the highest for the local and high-z star-forming galaxies. 
It is clearly seen that galaxies with younger starbursts (ma¬ 
genta symbols) have systematically higher H/3 luminosities, 
as compared to the galaxies with older starbursts (cyan sym¬ 
bols). 

For comparison, in Fig. 10a we show data for high-z 
star-forming galaxies by Cullen et al. (2014) and Steidel et 
al. (2014) (inside the box), Lya-emitting galaxies by Ha¬ 
gen et al. (2014) (stars), star-forming galaxies by Troncoso 
et al. (2014) (asterisks), star-forming galaxies with z ~ 4 
by Bouwens et al. (2012) (thick red solid line), star-forming 
galaxies with z ~ 5 - 6 by Salmon et al. (2015) (thick yel¬ 
low solid line), and galaxy candidates with z ~ 9 - 11 by 
Coe et al. (2013) and Oesch et al. (2014) (large open dark- 
grey triangles). High -z galaxies, including the highest red¬ 
shift galaxy candidates, nicely follow the trend for low -z 
galaxies with EW(H/3) ^100A (magenta symbols), but are 
offset to higher L(H/3) as compared to the rest of galax¬ 
ies from our sample with presumably older starbursts (cyan 
symbols). 

Reducing the H/3 luminosities (or, equivalently, SFRs) 
of galaxies from our sample to zero starburst age results 
in a much tighter relation in the redshift range 0 < z < 1 
(magenta and cyan symbols in Fig. 10b). This is because the 
age correction for galaxies with EW(H/3) ^ 100A is much 
smaller compared to that for the remaining galaxies of our 
sample. We note that the correction for the starburst age of 
the stellar mass is small because this parameter was obtained 
from the SED in the optical range, which is insensitive to 
age variations on the scale of several Myr. Therefore we did 
not introduce a correction for the stellar mass. 

It is seen in Fig. 10b that the distributions of low -z 
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Figure 10. The stellar mass - H/3 luminosity (or equivalently star-formation rate) relations for the entire sample. The extinction- and 
aperture-corrected H (3 luminosities and stellar masses are shown in panel (a), while H (3 luminosities are additionally reduced to zero 
starburst age in panel (b). Galaxies with the equivalent widths EW(H/3) <100A and ^100A are shown by the small cyan open circles and 
small magenta filled circles, respectively. Solid lines represent linear maximum-likelihood fits to all data. Constant specific star-formation 
rates (sSFR=SFR/M*) are shown by dotted lines and labelled by values in units of Gyr -1 . The location of star-forming galaxies with 
redshifts 2.0 < z < 2.6 by Steidel et al. (2014) and with z > 2 by Cullen et al. (2014) is shown by a grey rectangle. Lya emitting galaxies 
with redshifts 1.9 < z < 3.6 by Hagen et al. (2014) and star-forming galaxies with redshifts 2 : ~ 3.4 by Troncoso et al. (2014) are shown 
by dark-grey stars and asterisks, respectively. The relations for galaxies with 2 : ~ 4 by Bouwens et al. (2012) and for galaxies with 2 : ~ 
5 - 6 by Salmon et al. (2015) are indicated by thick solid red and yellow lines, respectively. The galaxy candidates with 2 : ~ 9 - 11 by 
Coe et al. (2013) and Oesch et al. (2014) are shown by large open dark-grey triangles. The H/3 luminosities of all high- 2 ; galaxies in (b) 
are not reduced to zero starburst age. 


and high- 2 ; star-forming galaxies in the Lo(H/3) - M*/M© 
plane are very similar, implying that the physical properties 
of these galaxies are similar as well. However, we note that 
H (3 luminosities for high -z galaxies are not corrected for the 
starburst age. 


6 SUMMARY 

We studied the relations between global parameters (abso¬ 
lute magnitudes, H f3 luminosities, star-format ion rates SFR, 


stellar masses, and gas-phase oxygen abundances) of 5182 
compact star-forming galaxies with high-excitation H II re¬ 
gions in the redshift range 0 < z < 1 selected from the SDSS 
DR7 and SDSS/BOSS DR10 surveys. These data were split 
into two samples, one is the T e sample of 3607 galaxies, 
where the temperature-sensitive [O ill] A4363 emission line 
was detected at the level better than lcr and the oxygen 
abundance was derived by the T e -method, and the second 
sample consisting of the remaining galaxies, where the oxy¬ 
gen abundances were derived by strong-line methods. The 
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data for local galaxies were combined with respective liter¬ 
ature data for high -z galaxies (with 2 < z < 3) to detect 
a possible redshift evolution of global parameters in star¬ 
forming galaxy. 

Our main results are as follows. 

1. Using the T e sample we analysed the oxygen abun¬ 
dances derived by the T e -method and various strong- 
emission-line (SEL) methods used for the chemical abun¬ 
dance determination. We adopted strong-line methods, 
which are consistent with the T e -method for high-excitation 
H II regions. 

2. We analysed various relations between absolute mag¬ 
nitudes, YL/3 luminosities (or equivalently star-formation 
rates SFR), stellar masses, and oxygen abundances for our 
sample of low -z compact star-forming galaxies in combina¬ 
tion with the respective data for high-z star-forming galax¬ 
ies. We find good agreement between relations for low -z and 
high -z galaxies indicating very weak redshift evolution of 
global parameters and weak dependence of metallicity on 
SFR, contrary to results in some other studies. This find¬ 
ing emphasizes the universal character of relations for con¬ 
sidered global parameters of compact star-forming galaxies 
with high-excitation H n regions at redshifts 0 < z < 3. 
Furthermore, galaxy candidates at z ~ 9 - 11 do not de¬ 
viate from low -z compact star-forming galaxies on the M* 
- SFR diagram implying that universal relations between 
global galaxy parameters may be valid for galaxies at higher 
redshifts. 


ACKNOWLEDGEMENTS 

Y.I.I., N.G.G. and K.J.F. are grateful to the staff of the 
Max Planck Institute for Radioastronomy (MPIfR) for their 
warm hospitality. Y.I.I. and N.G.G. acknowledge financial 
support by the MPIfR. Funding for the SDSS and SDSS- 
II has been provided by the Alfred P. Sloan Foundation, 
the Participating Institutions, the National Science Foun¬ 
dation, the U.S. Department of Energy, the National Aero¬ 
nautics and Space Administration, the Japanese Monbuk- 
agakusho, the Max Planck Society, and the Higher Edu¬ 
cation Funding Council for England. The SDSS Web Site 
is http://www.sdss.org/. The SDSS is managed by the As- 
trophysical Research Consortium for the Participating In¬ 
stitutions. The Participating Institutions are the American 
Museum of Natural History, Astrophysical Institute Pots¬ 
dam, University of Basel, University of Cambridge, Case 
Western Reserve University, University of Chicago, Drexel 
University, Fermilab, the Institute for Advanced Study, the 
Japan Participation Group, Johns Hopkins University, the 
Joint Institute for Nuclear Astrophysics, the Kavli Insti¬ 
tute for Particle Astrophysics and Cosmology, the Korean 
Scientist Group, the Chinese Academy of Sciences (LAM- 
OST), Los Alamos National Laboratory, the Max-Planck- 
Institute for Astronomy (MPIA), the Max-Planck-Institute 
for Astrophysics (MPA), New Mexico State University, Ohio 
State University, University of Pittsburgh, University of 
Portsmouth, Princeton University, the United States Naval 
Observatory, and the University of Washington. Funding for 
SDSS-III has been provided by the Alfred P. Sloan Foun¬ 
dation, the Participating Institutions, the National Science 
Foundation, and the U.S. Department of Energy Office of 


Science. The SDSS-III web site is http://www.sdss3.org/. 
SDSS-III is managed by the Astrophysical Research Con¬ 
sortium for the Participating Institutions of the SDSS- 
III Collaboration including the University of Arizona, the 
Brazilian Participation Group, Brookhaven National Lab¬ 
oratory, Carnegie Mellon University, University of Florida, 
the French Participation Group, the German Participation 
Group, Harvard University, the Instituto de Astrofisica de 
Canarias, the Michigan State/Notre Dame/JINA Partici¬ 
pation Group, Johns Hopkins University, Lawrence Berke¬ 
ley National Laboratory, Max Planck Institute for Astro¬ 
physics, Max Planck Institute for Extraterrestrial Physics, 
New Mexico State University, New York University, Ohio 
State University, Pennsylvania State University, University 
of Portsmouth, Princeton University, the Spanish Partic¬ 
ipation Group, University of Tokyo, University of Utah, 
Vanderbilt University, University of Virginia, University of 
Washington, and Yale University. 


REFERENCES 

Abazajian K. N. et ah, 2009, ApJS, 182, 543 
Planck collaboration: Ade P. A. R., Aghanim N., Armitage- 
Caplan C., Arnaud M. et ah, 2014, AVA, 571, A16 
Ahn C. P. et ah, 2014, ApJS, 211, 17 
Aller L. H., 1984, Physics of Thermal Gaseous Nebulae. 
Dordrecht: Reidel 

Alloin D., Coffin-Souffrin S., Joly M., Vigroux L., 1979, 
AVA, 78, 200 

Amorm R. O., Perez-Montero E., Vlchez J. M., 2010, ApJ, 
L128 

Andrews B. H., Martini P., 2013, ApJ, 765, 140 
Baldwin J. A., Phillips M. M., Terlevich R., 1981, PASP, 
93, 5 

Bouwens R. J., Illingworth G. D., Oesch P. A. et ah, 2012, 
ApJ, 754, 82 

Cardamone C. et ah, 2009, MNRAS, 399, 1199 
Cardelli J. A., Clayton G. C., Mathis J. S. 1989, ApJ, 345, 
245 

Chariot S., Longhetti M., 2001, MNRAS, 887 
Coe D., Zitrin A., Carrasco M. et ah, 2013, ApJ, 762, 32 
Cullen F., Cirasuolo M., McLure R. J., Dunlop J. S., Bowler 
R. A. A., 2014, MNRAS, 440, 2300 
Dawson K. S. et ah, 2013, AJ, 145, 10 
Dopita M. A., Sutherland R. S., Nicholls D. C., Kewley L. 
J., Fogt F. P. A, 2013, ApJS, 208, 10 
Fioc M., Rocca-Volmerange B., 1997, AVA, 326, 950 
Guseva N. G., Izotov Y. I., Thuan T. X., 2006, ApJ, 644, 
890 

Guseva N. G., Izotov Y. I., Papaderos P., Fricke K. J., 2007, 
AVA, 464, 885 

Guseva N. G., Papaderos P., Meyer H. T., Izotov Y. I., 
Fricke K. J., 2009, AVA, 505, 63 
Hagen A. et ah, 2014, ApJ, 786, 59 
Heckman T. M. et ah, 2005, ApJ, 619, L35 
Hoyos C., Koo D. C., Phillips A. C., Willmer C. N. A., 
Guhathakurta P., 2005, ApJ, 635, L21 
Hunt L. K. et ah, 2012, MNRAS, 427, 906 
Izotov Y. I., Thuan T. X., Lipovetsky V. A., 1994, ApJ, 
435, 647 


© 2012 RAS, MNRAS 000 , 1-13 



Universality of mass-metallicity relation 13 


Izotov Y. I., Thuan T. X., & Lipovetsky V. A., 1997, ApJS, 
108, 1 

Izotov Y. I., Stasinska G., Meynet G., Guseva N. G., Thuan 
T. X., 2006, A&A, 448, 955 

Izotov Y. I., Guseva N. G., Thuan T. X., 2011, ApJ, 728, 
161 

Izotov Y. I., Guseva N. G., Fricke K. J., Henkel C., 2014a, 
A&A, 561, A33 

Izotov Y. I., Guseva N. G., Fricke K. J., Kriigel E., Henkel 
C., 2014b, A&A, 570, A97 
Jaskot A. E., Oey M. S., 2014, ApJ, 791, L19 
Kakazu Y., Cowie L. L., Hu E. M., 2007, ApJ, 668, 853 
Kauffmann G., Heckman T. M., Tremonti C. et ah, 2003, 
MNRAS, 346, 1055 

Kennicutt R. C., Jr., 1998, Ann.Rev.Astron.Astrophys., 36, 
189 

Kewley L. J., Dopita M. A., 2002, ApJS, 142, 35 
Kewley L. J., Ellison S. L., 2008, ApJ, 681, 1183 
Kobulnicky H. A., Kewley L. J., 2004, ApJ, 617, 240 
Lara-Lopez M. A., Cepa J., Bongiovanni A. et ah, 2010, 
A&A, 521, L53 

Lee H. et ah, 2006, ApJ, 647, 970 
Leitherer C. et ah, 1999, ApJS, 123, 3 
Maier C., Lilly S. J., Ziegler B. L., Perez Montero E., Peng 
Y., Balestra I., 2014, ApJ, 792, 3 
Maiolino R. et ah, 2008, A&A, 488, 463 
Manucci F., Cresci G., Maiolino R., Marconi A., Gnerucci 
A., 2010, MNRAS, 408, 2115 

Nicholls D. C., Dopita M. A., Sutherland R., Kewley L. J., 
Palay E., 2013, ApJS, 207, 21 
Oesch P. A., Bouwens R. J., Illingworth G. D. et ah, 2014, 
ApJ, 786, 108 

Pagel B. E. J., Edmunds M. G., Blackwell D. E., Chun M. 
S., Smith G., 1979, MNRAS, 189, 95 
Palay E., Nahar S. N., Pradhan A. K., Eissner W., 2012, 
MNRAS, 423, L35 

Pettini M., Pagel B. E. J., 2004, MNRAS, 348, L59 
Pettini M., Shapley A. E., Steidel C. C. et ah, 2001, ApJ, 
554, 981 

Pilyugin L. S., 2001, A&A, 374, 412 
Pilyugin L. S., Thuan T. X., 2005, ApJ, 631, 231 
Pilyugin L. S., Vflchez J. M., Thuan T. X., 2010, ApJ, 720, 
1738 

Refsdal S., Stabell R., de Lange F. G., 1967, Mem. RAS, 
71, 143 

Salmon B., Papovich C., Finkelstein S. L. et ah, 2015, ApJ, 
799, 183 

Sanders R. L., Shapley A. E., Kriek M. et ah, 2015, ApJ, 
799, 138 

Smee S. A. et ah, 2013, AJ, 146, 32 
Stasinska G., 2005, A&A, 434, 507 
Stasinska G., 2006, A&A, 454, L127 
Steidel C. S. et ah, 2014, ApJ, 795, 165 
Tremonti C. et al., 2004, ApJ, 613, 898 
Troncoso P. et al., 2014, A&A, 563, A58 
York D. G. et al., 2000, AJ, 120, 1579 
Zahid H. J., Kewley L. J., Bresolin F., 2011, ApJ, 730, 137 
Zahid H. J., Bresolin F., Kewley L. J., Coil A. L., Dave R., 
2012, ApJ, 750, 120 

Zahid H. J., Geller M. J., Kewley L. J., Hwang H. S., 2013, 
ApJ, 771, L19 

Zahid H. J. et al., 2014a, ApJ, 791, 130 


Zahid H. J. et al., 2014b, ApJ, 792, 75 
Zhao Y., Gao Y., Gu Q., 2010, ApJ, 710, 663 

This paper has been typeset from a TgX/ DTgX file prepared 
by the author. 


© 2012 RAS, MNRAS 000 , 1-13 



